A Percolative Model of Soft Breakdown in Ultrathin Oxides 
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> ■ Abstract 

oo 

•£2 ' The degradation of ultrathin oxide layers in the presence of a stress voltage is modeled in 

O ' terms of two antagonist percolation processes taking place in a random resistor network. 

Q I The resistance and leakage current fluctuations are studied by Monte Carlo simulations 



for voltages below the breakdown threshold. An increase of excess noise together with a 



^ ■ noticeable non-Gaussian behavior is found in the pre-breakdown regime in agreement with 

■ experimental results. 
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1. Introduction and model 

Soft breakdown occurs in ultrathin oxide layers during constant voltage stress. This process 
manifests itself by a large increase of the stress-induced leakage currents (SILC) and by 
the occurrence of giant current fluctuations associated with resistance fluctuations. These 
fluctuations can be related to trapping and detrapping of electrons in the percolation cluster 
formed between the electrodes of the capacitor at the soft breakdown [1,2]. In this paper, 
we investigate the resistance in the pre-breakdown region and the associated fluctuations 
using two percolations which evolve in competition on a random resistor network [3], as 
briefly summarized below. We describe a thin fllm as a two-dimensional square-lattice 
network of resistors of resistance r„, laying on an insulating substrate at temperature 
To acting as a thermal reservoir. Geometrical effects arising from different shapes are 
neglected. Initially all the resistors are identical, r„ = tq. We take a square geometry 
N X N where N determines the linear sizes of the lattice and Nfot = SA?"^ is the total 
number of resistors. Electrical contacts, realized by perfectly conducting bars, are placed 
at the left and right hand sides of the network where a constant voltage V is applied. 
By increasing the applied voltage, defects are introduced by replacing single resistors with 
short circuits. Thus, the network degrades by undergoing an insulator-conductor transition 
when the defects gradually become so dense that flnally they form a continuous short circuit 
path between the contacts [3] . The temperature Tn of the n-th resistor is computed by [3] : 



where the term which adds to Tq in the r.h.s. of Eq. (1) accounts for Joule heating associ- 
ated with the current flowing in the resistor and the coupling with the nearest neighbours. 
Here, N^eig is the number of first neighbours around the n-th resistor with current z„, 
A the heat coupling parameter and 5 = 3/4 provides a uniform heating everywhere in 
the perfect network. We note that Eq. (1) assumes an instantaneous thermalisation be- 
tween neighbor resistors. Probabilities of creating and recovering a defect are taken as 
Wd = exp{—ED/KBTn) and Wn = exp{—Eji/KBTn) with and Eji characteristic ac- 
tivation energies, respectively, and Kb the Boltzmann constant. Monte Carlo simulations 
are carried out using the following procedure, (i) Starting from the perfect lattice, we 
calculate the total network resistance i?, the external current /, and the local currents by 
solving Kirchhoff's loop equations by the Gauss elimination method. The local tempera- 
tures are then calculated and used for the successive update of the network, (ii) The defects 
are generated with probability Wd- The local currents and the local temperatures are then 
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recalculated, (iii) The defects are recovered with probability Wr, and the total network 
resistance, the local currents, and the temperatures are finally calculated. This procedure 
is iterated from (ii). Each iteration step can be associated with an elementary time step. 
The iteration runs until either the defect percolation threshold is reached (i.e. the network 
resistance drops below 10~^ times to its initial value), or steady-state is achieved. In this 
latter case, the iteration is let to run long enough for a fluctuation analysis to be carried 
out. 

2. Results and discussion 

We compare the results of numerical simulations to the corresponding experimental ob- 
servations and validate the present approach by showing that it reproduces the key ex- 
perimental features. Since soft breakdown experiments report leakage current evolutions, 
noise spectra and probability density of fluctuations, we investigate the same quantities 
by numerical simulations. In the simulations we applied a constant voltage stress with 
the following realistic parameters N = 75, tq = 10^ fl, Tq = 300 K, A = 5 x 10^ °C/W, 
Ed = 0.19 eV, = 0.13 eV. Figure 1 reports the leakage current evolutions for stress 
voltages, respectively, of 0.1, 0.3, 0.5, 0.6, 0.7, 0.8, 0.85 and 0.9 V . We note that the values 
of the activation energies Ed and En control the speed of the evolution toward steady 
state and/or breakdown and as such should be calibrated from a quantitative compari- 
son with experiments. However, the qualitative features of the breakdown process we are 
presently interested in do not depend signiflcantly from their absolute values. We found 
that for voltages lower than 0.9 V no breakdown occurs but a dynamical balance takes 
place between defect generation and defect recovery. The steady state so achieved exhibits 
current fluctuations whose statistical and spectral features characterize the electrical qual- 
ity of the dielectric. However, for stress voltages equal or larger than 0.9 V after a few 
thousand iteration steps the leakage current grows more than a thousand fold which in- 
dicates a breakdown situation. We also observe that the breakdown does not occur as a 
result of a systematic current growth but more as a result of current bursts. In fact, a 
set of bursts generally occurs that lead to the breakdown of the insulator. Thus an inter- 
esting pre-breakdown region appears like in experiments [4,5]. From a detailed analysis 
of a set of simulations we also conclude that a threshold voltage exists below which no 
breakdown occurs. Figure 2 shows how the breakdown voltage depends on the resistance 
of the virgin dielectric, modeled as an initial perfect lattice with increasing value of the 
elemental resistor. This breakdown voltage is found to increase proportional to the square 
root of the perfect network resistance. Figure 3 shows the steady state current noise spec- 



trum calculated by FFT for stress voltages ranging from 0.1 to 0.85 V. Within numerical 
uncertainity, we detect Lorentzian spectra of the same corner frequency which indicates 
that the characteristic times of fluctuations are independent of stress voltages. Interest- 
ingly, the strong increase of the value of the plateaux at increasing stress voltages shows a 
super-quadratic increase of the noise with the applied stress. The current-noise spectra in 
Fig. 3 are found to be in qualitative agreement with experiments [4] . Figure 4 reports the 
normalized variance of resistance fluctuations as a function of the applied voltage. Starting 
from the intrinsic value of the network [6], the variance is found to increase signiflcantly 
as a net effect of the stressing voltage which is ultimately responsible of the breakdown 
at about 0.9 V. Recent observations of SILC measurements reported non-Gaussian cur- 
rent fluctuations in the soft breakdown region of ultra-thin dielectrics [1]. Therefore, we 
computed the distribution function of the current fluctuations P{I) for stress voltages of 
0.5 V and 0.85 V and plotted them in Fig. 5 as a function of / — < / >, where < / > is 
the average current value. On this linear-log plot the corresponding Gaussian distribution 
would be a downward parabola. The distribution for an applied voltage near breakdown 
of 0.85 V clearly exhibits a non-Gaussian tail at high currents over the average due to the 
large upward bursts which are not balanced by downward bursts in the current evolution. 
Since a further increase of the stress voltage causes the fllm to break down, the emerg- 
ing non-Gaussian distribution can be thought of as a precursor of failure as observed in 
experiments [1]. 

3. Conclusions 

We have developed a percolative approach to study the soft breakdown and the asso- 
ciated current fluctuations of ultrathin insulating films. The major results we found are 
summarized as follows, (i) The current evolution exhibits a pre-breakdown region, (ii) 
A threshold voltage depending on the resistance of the otherwise perfect film separates 
breakdown and steady state conditions, (iii) The spectrum of the current fluctuations 
is basically of Lorentzian type with a correlation time which is practically independent 
of the stress voltage, (iv) The variance of resistance fluctuations exhibits a giant super- 
quadratic enhancement at increasing the stress voltage, (v) At high enough voltage stress 
the current fluctuations exhibit a non-Gaussian behavior. These features, particularly 
the super-quadratic and the non-Gaussian behaviors of the current noise at high stress, 
are of relevant interest as failure precursors. These results are in satisfactory qualitative 
agreement with experimental findings. 
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Figure Captions 



Figure 1 - Leakage current evolutions in a degrading dielectric kept under constant voltage 
for stress voltages (from bottom to top) of: 0.1, 0.3, 0.5, 0.6, 0.7, 0.8, 0.85 and 0.9 V. 

Figure 2 - Breakdown voltage Vi, as a function of the perfect network resistance Rper- The 
graph shows oc yjR^r- 

Figure 3 - Power spectral density of current fluctuations of the steady-state dielectric film 
kept under constant voltage of (from bottom to top): 0.1, 0.3, 0.5, 0.6, 0.7 and 0.85 V . 

Figure 4 - Normalized variance of the network resistance fluctuations as a function of the 
applied voltage. 

Figure 5 - Gaussianity check of current fluctuations for a stress voltage of 0.5 V (open 
squares) and of 0.85 V (flUed circles). 
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